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ABSTRACT 

RX J0720.4— 3125 is an unidentified bright soft X-ray source which shows pulsations 
at a 8.39 s period and has a thermal spectrum. We present deep B and R band images 
of its X-ray localization. We find one possible counterpart in the X-ray error box, with 
magnitudes B = 26.6 ±0.2 and R = 26.9 ±0.3. The very high X-ray to optical flux ratio 
confirms that this object is an isolated neutron star. We discuss possible models and 
conclude that only two are consistent with the data and at the same time are able to 
draw from a large enough population to make finding one nearby likely. In our opinion 
the second criterion provides a stringent constraint but appears to have been ignored so 
far. The first model, suggested earlier, is that RX J0720.4— 3125 is a weakly magnetized 
neutron star accreting from the interstellar medium. The second is that it is a relatively 
young, highly magnetized neutron star, a "magnetar", which is kept hot by magnetic 
field decay. 

Subject headings: stars: individual (RX J0720.4— 3125) — stars: neutron — stars: mag- 
netic fields — X-rays: stars 
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1. Introduction 



2. RX J0720.4 3125 



The population of defunct radio pulsars far exceeds 
that of active ones. It is believed that the Galaxy has 
about 2 10 5 radio pulsars. The neutron-star birth rate 
is estimated to be between one per 30 yr to one per 
100 yr. Assuming a constant pulsar production rate 
and an age of the disk of 10 10 yr, one infers a Galactic 
neutron star population of ~ 2 10 s - three orders of 
magnitude larger than that of the active radio pulsar 
population. 

It is not easy to detect old neutron stars. While 
the nearest few intermediate-age pulsars can be iden- 
tified by their cooling radiation, which peaks in the 
soft X-ray/EUV band, the defunct pulsars will have 
become too cool to be observable. A small fraction of 
them, however, may be in a position to accrete mat- 
ter from the interstellar medium. These will then get 
reheated and reappear in the X-ray sky. 

Quite independent of this discussion there has been 
a growing recognition of a population of highly mag- 
netized neutron stars. The circumstantial evidence 
for this class comes from studies of soft gamma-ray 
repeaters (SGRs) and long-period pulsars in super- 
nova remnants ( Vasisht & Gotthelf 1997| ) . Thompson 
& Duncan (1995) have introduced the term "mag- 
netars" for neutron stars with field strengths signif- 
icantly larger than 10 12 G, the typical field strength 
inferred for radio and X-ray pulsars. The birthrate of 
SGRs has been estimated to be roughly 10% that of 
the ordinar y puls ars (Kulkarni & Frail 1993 ; Kouve- 
liotou et al. 



1994| ) 



The relevance of magnetars to the discussion at 
hand is as follows. Unlike the situation for ordinary 
neutron stars, magnetic field decay is expected to be 
significant in highly magnetized neutron stars. This 
dec ay cou ld reheat the magnetar (Thompson & Dun- 



can [1996 ) , making it hotter than an ordinary neutron 
star, and thus brighter in soft X rays. 

Two of the best candidates for this general class of 
neutron stars have emerged from the ROSAT mission: 
RX J185635-3754 (Walter, Wolk, & Neuh auser |l99^ ) 
and RX J0720.4-3125 (Haberl et al. |l997j ). Both are 
bright ROSAT objects with very soft X-ray spectra. 
Walter & Matthews ( 1997| ) have provided compelling 
evidence for the identification of a faint blue optical 
counterpart of RX J185635— 3754. In this Letter, we 
present deep B and R observations of the localization 
of RX J0720.4-3125. 



The basic X-ray properties are described by Haberl 
et al. ( 1997 ). A brief summary now follows. With a 
PSPC count rate of 1.67 s _1 the source is among the 
brighter in the ROSAT All-Sky Survey. The PSPC 
spectrum is fitted very well by a black body with 
fcT ff = 79±4eV absorbed by a column density, iV# = 
(1.3 ± 0.3) 10 20 cm- 2 . The observed flux in the 0.1- 



2.4 keV ROSAT band is 1.15 



+0.3 
-0.14 



io- 



erg cm s 



Using PI MMS (Portable, Interactive Multi-Mission 
Simulator; http://legacy.gsfc.nasa. gov ), we estimate, 
for the black-body model mentioned above, an un- 
absorbed flux of 1.8 10" 11 erg cm" 2 s" 1 in the range 
0.010-2.4 keV. 

Perhaps the most remarkable fact about the source 
is neither its high count rate nor its soft X-ray spec- 
trum, but the presence of pulsations with a 8.39 s pe- 
riod. The period is very stable: from observations 
done over a three year period Haberl et al. find a 
period derivative P = (-2.6 ± 3.4) lO^ss" 1 . This 
argues for the period to be the rotation period of the 
star. The shortness of the period rules out all models 
save those involving a neutron star or a white dwarf. 

Haberl et al. found no optical counterpart within 
the ROSAT HRI error circle and noted a X-ray to 
optical flux of > 500. Such a high value is barely 
compatible with a model involving an accreting low 
mass X-ray binary and certainly excludes any model 
based on an accreting white dwarf binary. We note 
that an isolated hot white dwarf (planetary nebula 
nucleus) would be incompatible with the high inferred 
black-body temperature. 

3. Optical Observations 

We imaged the field containing RX J0720.4-3125 
with the Low Resolution Imaging Spectrograph (Oke 
et al. 1995), mounted at the Cassegrain focus of the 



Keck II telescope, on the nights of November 28 and 
29, 1997 (UT). The detector was a 2048 x 2048 pixel 
Loral CCD which was read out with two amplifiers. 
The skies appeared to be clear on the second night 
but the first night was plagued by cirrus. A log of the 
observations is given in Table @. 

During the second night, four flux-standard fields 
from Landolt Ql992| ) were observed (defocusing the 
telescope to avoid saturation of pixels). We took care 
to place our targets and most of the standards on the 
same side of the CCD chip, thereby ensuring that the 
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same amplifier was used. From the standard stars, 
which had a range in B — R from —0.4 to 1.7, we 
determined the color and extinction coefficients. The 
root-mean-square residual around the calibration is 
about 0.02 mag in both B and R. 

The reduction was done using MIDAS. Individual 
frames were bias-subtracted and then flat-fielded us- 
ing twilight flats. For the RX J0720.4-3125 images, 
a fixed region around the target's position was ex- 
tracted and cleaned of cosmic rays, and these images 
were added together by band and by night. The result 
is shown in Figure [l] (Plate) . Overdrawn are the two 
X-ray position estimates from Haberl et al. (1997), us- 
ing astrometry relative to the USNO-A1.0 catalogue 
(Monet et al. |1996| ); see also Table |. Note that these 
are close to, but not exactly at, the position shown in 
Fig. 7 of Haberl et al. We do not know the reason for 
this discrepancy. 



As can be seen in the figure, only one source (here- 
after, "X") is found within the union (and intersec- 
tion) of the two HRI error circles. It is well detected in 
the B image but barely detected in the R band. Scat- 
tered light from some nearby stars result in a strong 
gradient in the local sky. This made application of 
the standard photometric packages difficult. We pro- 
ceeded by using a simplified point-spread function fit- 
ting, as follows. First, we used aperture photometry 
to measure the magnitudes for four relatively isolated, 
brighter "secondary" stars (22, 24, A, B; see Figure [y 
and Table ||). Next, we extracted the stellar images 
of these four stars, as well as those of the candidate 
X and of star C, within a 21 x 21 pixel (4'.'5 x 4"5) 
region centered on each source, and fitted these to 
a 2-dimensional Gaussian on top of a plane with an 



Table 1 
Log of Observations 



Table 2 

Photometry and Astrometry of the 
Candidate Counterpart of RX J0720. 4-3125 
and Selected Stars in the Field 



Field 


Filter 


UT 


^cxp 


sec z 


Star 


B 


R 


OJ2000 


£12000 


RX J0720.4-3125 


R 


12:44 


2 x 300 


1.60 


X 


26.6(2) 


26.9(3) 


7 20 24.92 


-31 25 50.9 




B 


12:58 


4 x 600 


1.60 


A 


20.52 


19.42 


7 20 24.40 


-31 25 52.3 




R 


13:47 


3 x 300 


1.63 


B 


20.79 


19.62 


7 20 25.15 


-31 26 00.3 












C 


21.23 


18.85 


7 20 25.76 


-31 25 55.6 


RX J0720.4-3125 


B 


12:27 


3 x 900 


1.60 


22 


20.37 


19.22 


7 20 24.18 


-31 26 28.4 




R 


13:18 


7 x 300 


1.62 


24 


21.49 


20.13 


7 20 22.88 


-31 26 04.6 


PC 2336+004 


B 


08:28 


20 


1.50 












(AM) 


R 


08:32 


5 


1.52 












SA 95 


R 


08:38 


5 


1.10 












(41,42) 


B 


08:45 


20 


1.10 


Note. — The uncertainty in the photometry for 


PG 1047+003 


R 


14:08 


5 


1.33 


star X 


is indicated by the number in 


brackets. For 


(PG,A,B,C) 


B 


14:10 


20 


1.33 


the other stars, 


it is <; 0.02 mag in both B and R. 


Ruben 149 


B 


15:43 


10 


1.35 


Comparing our measured B — R to the 


main sequence 


(Ru,B,C,E,F,G) 


R 


15:45 


2x2 


1.36 


colors (Johnson 


1966) we 


conclude that the spectral 



Note.— The UT date is 1997 Nov. 28 for the first 
three entries, 1997 Nov. 29 for all others. All exposure 
times are in seconds. The standard fields are from 



types of A and B are both ~ G2 - quite consistent 
wit h wha t is indicated by the spectra of Haberl et 
al. 



(1997) 



Landolt (1992); below the name the stars actually 
used for the photometric calibration are listed. The 
seeing (as measured in the R-band) on the first night 
was 1" and improved to 0"8 on the second night. 



Haberl et al. classified star C as M based 
on its spectrum - from both color and spectrum, we 
believe it is a little earlier, ~ K7. The V magnitudes 
inferred using these spectral types are consistent with 
the range quoted by Haberl et al. The uncertainty in 
the astrometry is dominated by the systematic uncer- 
tainty in the USNO-Al.O catalogue. 
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20 m 28 s 20 m 24 s 20 m 25 s 20 m 24 s 



«2000" 7h 

Fig. 1. — Summed B (top) and R (bottom) images of the localization of RX J0720.4— 3125. The left-hand panels 
show an overview. The two circles represent the two HRI X-ray positions derived by Haberl et al. (1997). Stars 
mentioned in the text are labeled below their image in the top panel. For each band, the region around the localiza- 
tion is shown enlarged on the right-hand side. For these enlargements, the local sky was subtracted (determined by 
fitting a two-dimensional, second-order polynomial to the image with stars masked out). In the lower panels, the 
image has furthermore been convolved with a Gaussian with a width approximately matching the seeing (FWHM 
of four pixels - Of! 85). 
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arbitrary tilt. From the fit, we determined the av- 
erage FWHM for the four secondary stars. We then 
refitted all objects keeping the FWHM fixed at the 
average, and used the amplitudes of the Gaussians to 
determine relative magnitudes. Finally, the magni- 
tudes of stars X and C were measured relative to the 
secondary stars which in turn were calibrated using 
the solution found from the secondary stars. 

In summary, our observations have identified only 
one source, star X, within the HRI localization re- 
gion. This source is a plausible optical counterpart of 
RX J0720.4-3125. Its magnitudes are B = 26.6 ±0.2 
and R = 26.9 ± 0.3 (Table |). It is moderately blue, 
B-R = -0.3 ±0.4. 

4. The Nature of RX J0720.4 3125 

A straightforward conclusion we can draw from our 
observations is that the ratio of the X-ray to the opti- 
cal fluxQ, fx/fv <; 210 5 . The approximate equal- 
ity is applicable if X is the optical counterpart of 
RX J0720.4-3125. The high value for the ratio rules 
out the possibility that RX J0720.4— 3125 is an accret- 
ing neutron star binary, leaving only a single neutron 
star as a viable option. The X-ray spectrum is highly 
suggestive of thermal emission from a hot surface. If 
so, fx constrains the distance to the source to be less 
than 450i?6pc; here Rq is the radius of the neutron 
star in units of 10 6 cm. 

RX J0720.4— 3125 must be nearby since the in- 
ferred column density is rather small, 1.3 10 20 cm~ 2 . 



As pointed out by Haberl et al. (1997), the source is in 
the general direction of the o pen cluster Collindcr 140 
(Claria fc R osenzweig 197S). FitzGerald, Harris & 
Miller ( |l980|) carried out extensive spectroscopic ob- 
servations and derive Eb-v — 0.04, a distance of 
410±30pc, and an age of 20±6Myr. Using the usual 
relation between dust and gas, the inferred total col- 
umn density is 2.8 10 20 cm- 2 . RX J0720.4-3125 is 
within the confine s of t he cluster , and thus, follow- 
ing Haberl et al. ( 1997| ), we place an upper limit of 
410 pc on the distance to RX J0720.4-3125. 

Accepting the conclusion that RX J0720.4— 3125 is 
an isolated and nearby neutron star, we now turn to 
the question of why its surface is hot. 



computing the X-ray to optical flux it is traditional to 

compute the optical flux using the following relati on: J y 
10 -0.4(V+13.42) ergcm -2 f 

assumed V ~ B for source X. 



see Maccacaro et al. (1988). We 



4.1. A Young Neutron Star 

In this model, RX J0720.4— 3125 is a nearby youn g 
neutron star. According to Umeda et al. ( |1993| ), 
an isolated neutron star even with internal frictional 
heating will cool down to kT = 80 eV at an age 
Tcooi = 4 10 5 yr. For a given neutron star birthrate 
-B ns , the expected number of neutron stars with an 
age less than r coo i and within a distance projected 
on the Galactic plane less than /9ii m is B na T coo \T:pf iai . 
We make the assumption that most neutron stars are 
born as ordinary pulsars. Lyne et al. ( |1998| ) estimate 
the local birth rate of pulsars with a 400-MHz lu- 
minosity above 1 mJykpc 2 as 2.8 ± 1.7Myr _1 kpc -2 . 
Taking into account that many pulsars will not be 
beamed towards us, they derive a total birth rate 
of 10 ± 6 Myr -1 kpc -2 . Lyne et al. claim that the 
luminosity function flattens below 20mJykpc 2 . If 
so, B ua is well constrained. Hence, for t coo i = 
0.4 Myr, pn m = 0.4 kpc, we expect a total of 1.2- 
3.2 young nearby pulsars showing strong thermal 
emission. Geminga, at 160 pc, is perhaps the near- 
est example of a mi ddle-aged cooling neutr on star 
(Caraveo et al. |l996| ; Halpern & Wang |1997| ), while 
PSR B0656+14 (Finley, Ogelman, & Kizilo glu |1992[ ) 
and PSR B1055-52 (Ogelman & Finley |l993j ), both at 
~ 500 pc, are examples of neutron stars whose radio 
emission is beamed towards us. 

A problem for the young-pulsar model, however, 
is the very long spin period of 8.39 s. There is lit- 
tle evidence that normal pulsars are born with such 
long periods, and it would have taken too long to spin 
down from an initially short period. Thus, we con- 
sider this model to be unlikely. 

4.2. A Neutron Star Accreting from the ISM 

In this picture, already extensively discussed in 
the literatur e (Ko nenkov & Popov |l997| ; Wang [1997 ; 
Haberl et al. |1997| ), it is assumed that the neutron star 
is heated by matter being accreted from the ISM. An 
advantage compared to the young pulsar model is that 
one can draw upon the entire old neutron star popula- 
tion. The principal uncertainty is whether all the con- 
ditions for accretion are met. These include: (1) the 
unknown and small filling factor of moderately dense 
interstellar medium; (2) the small fraction of neutron 
stars with sufficiently low velocity (<; 10%; Cordes 
& Chernoff 1997); and (3) the combination of dipole 
field strength and spin period needed to accrete mat- 
ter. Accretion can proceed only if the Alfven radius 
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(the radius at which accretion flow pressure balances 
the magnetospheric pressure) is inside the corotation 
radius. For RX J0720.4-3125, Wang Ql997j ) infers 
B < 10 10 G. Thus, if the field originally was similar 
to that observed in present-day young neutron stars, 
it must have decayed substantially. There is little 
evidence for any decay in ordinary pulsars. Indeed, 
Kulkarni & Anderson ( 1996| ) argue that the existence 
of long period pulsars in clusters is incompatible with 
a field decay by more than a factor of 10 over a Hubble 
time. For this reason we neglect this model. 



Haberl et al. (1997) suggested an interesting pos- 
sibility that might alleviate points (2) and (3) men- 
tioned above, namely that RX J0720.4— 3125 is an 
evolved member of the group of so -calle d anomalous 
X-ray pulsars (Mereghetti & Stella 1995 ). These sys- 
tems have similar pulse periods, and are found close 
to the plane, often associated with young SNRs. If 
these systems are formed in common-envelope evo- 
lution of high-mass X-ray binaries, as suggested by 
van Paradijs et al. ( 1995 ), their space velocities are 
likely low, and their magnetic fields may be reduced 
due to accretion. For a birth rate of anomalous pul- 
sars of about one per 3 10 3 yr (van Paradijs et al. 
1995 ), there should be ~ 1.6 10 3 descendants within 
P\im < 0.4 kpc. If all of these have suitably low veloc- 
ity, the fraction visible as X-ray sources will depend 
only on the filling factor / of dense (;> lcm~ 3 ) in- 
terstellar medium. Even for an unrealistic / ~ 0.01 
we expect to see many accreting systems. Thus, this 
is a plausible model, subject only to the assumed 
birthrate of anomalous pulsars and whether they in- 
deed evolve to low velocity, low magnetic field, single 
neutron stars. 

4.3. A Magnetar 

The long period of 8.39 s is expected if the ob- 
ject is not a normal young pulsar, but a magnetar. 



This possibility is briefly mentioned by Wang (1997) 
and c onsidered in more detail by Heyl & Hernquist 
( 1998 ). In the framework of the rotating magnet 



(vacuum) model for pulsars, one has PP oc B^; here 
P is the rotation period and B& the strength of the 
dipolc component of the magnetic field. Assuming 
that the neutron star was born with a period much 
shorter than the current period we find P = P/t = 
3 10 _13 r 6 " 1 s s _1 where the age of the magnetar is 
t = 10 6 t 6 yr. This leads to B d ~ 5t~ 1/2 10 13 G. 
Earlier we had stated that soft gamma-ray re- 



peaters have been argued to be young magnetars 
and that their birthrate is about 10% of that of 
the ordinary pulsars. With one magnetar for ev- 
ery ten ordinary pulsars, the probability of finding a 
middle-aged magnetar is small but not negligible. Ob- 
servationally, if we associate RX J0720.4— 3125 and 
RX J185635— 3754 with magnetars, the local number 
density of magnetars would be comparable to that 
of middle-aged ordinary pulsars. This would indi- 
cate that magnetars are brighter and/or longer -lived 
as soft X-ray sources. Heyl and Hernquist ( 1998 ) 
have studied the thermal evolution in the presence 
of magnetic fields in detail. Curiously, their cool- 
ing timescales appear insensitive to the field strength 
when kT e g <J 100 eV. However, their models ne- 
glect the influence of magnetic-field decay. For these 
high fiel d stre ngths, decay is expected (Thompson & 

Reisenegger |!992f) 



Duncan 1996 



also Goldreich 
We believe this field decay is what sets magnetars 
apart from and could make them brighter than ordi- 
nary neutron stars. Therefore, we would argue that 
RX J0720. 4-3125 could well be a magnetar. 

5. The nature of star X 

We now discuss the nature of star X. There are 
two possibilities. The first is that X is unrelated to 
RX J0720.4-3125. If so, X has to be a distant star, 
with spectral type no later than F2 (B — R < 0.72; 
3(7 limit). Combined with the apparent B magnitude 
of 26.6, it could only be a white dwarf (a main se- 
quence star would be well outside even the extended 
halo of our Gala xy) . A white dwarf would have 
My <; 13.2 (Allen |197§| ), and be at a distance ^4kpc, 
;>0.6 kpc above the plane. It would have to be a disk 
white dwarf, since white dwarfs in the halo would 
be too cool. The scale height of disk white dwarfs is 
275±50pc, and the local density 0.6 10~ 3 pc -3 (Boyle 
1991 ). At th e latitude of RX J0720.4-3125, -7?8, we 
would thus expect about one white dwarf per 4 square 
arcmin. Hence, the probability of finding one within 
the 80 square arcsec HRI error box is ~ 6 10 -3 . We 
conclude that it is possible, though not very likely, 
that X is a distant white dwarf. A better color de- 
termination could settle this issue; if the object really 
has B — R < 0, the putative white dwarf would have 
to be much brighter and thus be well outside the disk. 

This brings us to our second hypothesis, that star 
X is the optical counterpart of RX J0720. 4-3125. 
The observed B and R magnitudes correspond to 
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fluxes of 92 ± 14 and 53 ± 17nJy, respectively. Here, 
we used the effective wavelengths and absolute cali- 
bration given by Bessel ( |1992| ): X oS (B,R) = 0.436, 
0.638 jim, F U (B,R = 0) = 4.00, 3.06 kJy. We can 
compare these fluxes with those expected from the 
ROSAT data, by extrapolating the best-fit, kT g — 
80 eV, black-body spectrum (see Sect. ||). Using the 
Rayleigh-Jeans formula, f v = / c) 2 '{R/ 'd) 2 , 

and (R/d) 2 ~ fx/&T^ s - where fx is the unabsorbed 
flux in the ROSAT band (Sect, g) - the predicted B 
and R-band fluxes are 18 and 8 nJy, respectively. The 
observed optical flux is a factor of five larger. A sim- 
ilar excess is seen in RX J185635— 3754 (Walter & 



Matthews |1997| ). 

There are several possible causes for the apparent 
excess. Perhaps the simplest is that the surface does 
not have a uniform temperature; a two-component 
model, one component with kT\ — 80 eV and the 
other with kT 2 <. 80 eV but occupying 4Ti /T 2 times 
more area, would satisfactorily account for both the 
optical and the X-ray measurements. An alterna- 
tive is that the excess is only apparent; the emergent 
X-ray spectrum depends on the composition in the 
neutron-star atmosphere and the effective tempera- 
ture may thus well be different from that inferred by 
fitting a black-body spectrum (Pavlov et al. 1996). Fi- 
nally, the optical excess could be due to non-thermal 
emission, similar to what has been inferred for PSR 



B0 656+ 14 (Pavlov, W elty, fc Cordova |1997| ; Kurt et 
al. 1997; Shearer et al. 1997) and Geminga (Bignami 
et al. |1996| ). 

The measurement of P would clearly be the most 



ele4ant way to discriminate between the two mod- 



els tor KA J0720. 4— '612b. For the magnctar, one 
would expect steady spin-down, at a few 10 - ss -1 
(Sect. p~5|). For the accreting old neutron star model, 
Wang (1997) estimates a spin-down rate <^ 10~ 16 s s _1 . 



This ignores the torques exerted by the accreting mat- 
ter, however, which might lead to either spin-up or 
spin- down, at a much larger rate (Lipunov & Popov 
1995 ). Other future observations (improved position 
and better spectrum, AXAF; accurate optical and UV 
photometry, HST; radio emission, VLA) will be of 
great interest no matter what the object turns out to 
be. 

While revising this manuscript, we became aware 
of a preprint by Motch & Haberl ( 1998 ) reporting op- 
tical observations of RX J0720.4— 3125. They propose 
two candidates, XI and X2. XI is our star X. X2 is 
about 3 arcseconds NW of XI (see Figure [0) and is 



nominally outside the HRI error circles based on our 
astrometry. It is easily detected in R but fainter than 
XI in B. Given the abundance of such faint red objects 
it is our opinion that X2 is most likely a background 
object. 
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